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We investigate quantum phase transitions among the spin-gap phases and the magnetically 
ordered phases in a two-dimensional frustrated antiferromagnetic spin system, which interpolates 
several important models such as the orthogonal-dimer model as well as the model on the 1/5- 
depleted square lattice. By computing the ground state energy, the staggered susceptibility 
and the spin gap by means of the series expansion method, we determine the ground-state phase 
diagram and discuss the role of geometrical frustration. In particular, it is found that a RVB-type 
spin-gap phase proposed recently for the orthogonal-dimer system is adiabatically connected to 
the plaquette phase known for the 1/5-depleted square-lattice model. 
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§1. Introduction 

Geometrical frustration plays an important role for 
the spin-gap formation in a certain class of antiferro- 
magnetic quantum spin systems. Two typical spin- 
gap compounds found in two dimension (2D) are the 
transition metal oxides, SrCu2(B03)2^'' and CaV409.^-' 
In the former compound, the orthogonal-dimer struc- 
ture of the Cu^+ spins stabiUzes the dimer phase with 
the spin gap, which is properly described by the 2D 
Heisenberg model on the square lattice with diagonal in- 
teractions (Shastry-Sutherland model), '^^ as claimed by 
Miyahara and Ueda.^^ For the model, it was pointed 
out that there should exist an intermediate phase be- 
tween the dimer and the antiferromagnetic phases, 
implying that strong frustration due to the competing 
interactions is quite important for understanding this 
model. On the other hand, the latter spin-gap com- 
pound CaV409 may be described by the 2D Heisen- 
berg model on the 1/5-depleted square lattice^-' or the 
meta-plaquette model. ^' The spin-gap phase in this 
compound results from the plaquette-singlet formation, 
which is further stabilized by the frustration effect due 
to the next-nearest- neighbor interactions. 

An interesting feature common to the above systems is 
that geometrical frustration due to the competing inter- 
actions plays an essential role to stabilize the disordered 
ground state. It is thus important to deal with these sys- 
tems in a unified framework to clarify the role of frustra- 
tion systematically. In particular, since a RVB-type spin- 
gap phase^^ proposed for the Shastry-Sutherland model 
may be caused by geometrical frustration, it is desired 
to further clarify the nature of this spin-gap phase. This 
is also important from the experimental point of view, 
because the compound SrCu2 (603)2 is expected to be 
located around the phase boundary between the dimer 
and RVB-type phases. 

Motivated by the above hot topics, we study the com- 



petition among the spin-gap phases as well as the mag- 
netically ordered phases in a 2D frustrated quantum spin 
system on the square lattice with some diagonal cou- 
plings. The basic idea in this paper is the adiabatic con- 
tinuation, which allows us to relate apparently different 
states which belong to the same phase at the fixed point. 
By using this model, we systematically describe quan- 
tum phase transitions for the orthogonal-dimer model 
and the Heisenberg model on the 1/5-depleted square 
lattice, where the ground-state phase diagram was al- 
ready discussed in detail. In particular, we show 
that a RVB-type spin-gap phase naturally emerges in the 
Shastry-Sutherland model, by observing the adiabatic 
evolution of the spin-gap state when the competing in- 
teractions are varied. 

This paper is organized as follows. In §2, we introduce 
a 2D frustrated spin model, which allows us to treat the 
above spin-gap compounds systematically. After briefly 
summarizing the series expansion method in §3, we inves- 
tigate the first- and second-order quantum phase tran- 
sitions among the spin-gap phases and the magnetically 
ordered phases in §4. We there determine the phase di- 
agram and clarify the role of geometrical frustration in 
this class of quantum spin systems. A brief summary is 
given in §5. 

§2. Model Hamiltonian 

We consider the Heisenberg model on the square lat- 
tice with some diagonal couplings, which is described by 
the Hamiltonian 



Si • Sj + J' ^ Si • Sj , 



(2.1) 



diagonal 



square 



where Si is the s = 1/2 spin operator at the z-th site, and 
J, J'l and J2 represent the antiferromagnetic exchange 
couplings. By putting J on the orthogonal dimer bonds 
(dashed lines), we then distribute the other couplings J( 
and J2 in two ways as shown in Figs. 1(a) and 1(b). In 
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the case (a), there is a diagonal bond J in each plaquette 

formed by the coupling J[ (bold line), while in the case 
(b) there is no diagonal bond in shaded plaquettcs. What 
is interesting in this generalization is that we can discuss 
several different models in the same framework. Namely, 
both models of (a) and (b) are related to each other 
via the Shastry-Sutherland model with J{ = Jj, and 
moreover they are reduced to the well-known models in 
the limiting case of J2 = 0: the orthogonal-dimcr chain 
for (a) and the spin system with a windmill structure for 
(b). Note that the latter is topologically equivalent to the 
Heisenberg model on the 1/5-depleted square lattice for 
the compound CaV409. Therefore, by using the above 
models we can systematically study how the geometrical 
frustration in the 2D model affects the generation of the 
spin gap. For later convenience, we introduce the ratios 
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Fig. 1. Frustrated spin models on the square lattice with some 
diagonal bonds. The bold, thin and dashed lines represent the 
coupling constants , Jj and J. The initial configuration of the 
isolated dimers (plaquettes) for the series expansion is indicated 
by the dashed lines (shaded squares). 



§3. Series Expansion Method 

In order to study quantum phase transitions, we make 
use of the series expansion method. ^"^^ This method is 
powerful especially for quantum spin systems with frus- 
tration, where quantum Monte Carlo simulations may 
suffer from sign problems. In the present context, we 
wish to mention that this method was successfully ap- 
plied to the meta-plaquette system,^' the orthogonal 
dimer system,^' ^' '^'^^ etc. To determine the phase dia- 
gram, we consider two initial configurations in the series 
expansion. The first one is a plaquette-singlet config- 
uration composed of the shaded squares in Figs. 1(a) 
and 1(b). In this case, we parameterize the antiferro- 
magnetic couplings labeled by the bold, thin and dashed 
lines as J{, J2 = AfcJ( and J = AJ{/q!, where we have 
introduced an auxiliary series-expansion parameter A in 
addition to the dimensionless couplings fc and a. The 
original model we are interested in is recovered by setting 



A = 1. Note that a RVB-type spin-gap phase, which 

has the spatially-extended disordered ground state, may 
be reached in this approach. The second one is a dimer 
configuration, by which we can describe the dimer phase. 
A convenient parameterization in the dimer expansion is 
J(= 1), J{ = \a and J'^ = \ka. In this way, these two 
different approaches enable us to perform the series ex- 
pansion starting from isolated dimer and plaquette sin- 
glets (A = 0). 

To carry out the series expansion explicitly, the origi- 
nal Hamiltonian eq.(2.1) is divided into two parts. 



H = Ho + XHi 

= X ^ S j • Sj -|- A ^ ^ Fjj S j • Sj 



(3.1) 



where x = J (J{) and Fy = a or ka {k or a~^) for 
the dimer (plaquette) expansion. The first term Hq is 
the unperturbed Hamiltonian which stabilizes the iso- 
lated dimer or plaquette singlets, whereas the perturbed 
Hamiltonian Hi labeled by A connects these isolated 
dimer or plaquette singlets, yielding a 2D frustrated spin 
system. We compute the staggered susceptibility xaf, 
the spin-triplet excitation energy E{q), and the ground 
state energy Eg as a power series in A. To estimate the 
susceptibility, we introduce the Zeeman term 

H' = h J2st-J2st , (3.2) 
jeA ieB 

where h is the staggered magnetic field and A{B) denotes 
one of the two sublattices. 

We determine the quantum phase transition point in 
three ways: (i) comparison of the ground state energy, 
(ii) divergence of the staggered susceptibility, (iii) van- 
ishing point of the spin gap. The latter two schemes may 
be applied to the second-order phase transition, while the 
first one is more effective to study the first-order phase 
transition. In order to obtain the accurate phase bound- 
ary, we further use the Dlog Fade approximants^^^' for 
the calculated results up to the finite order in A. We also 
make use of the biased Fade approximants,^^) in which 
the critical exponents for our 2D quantum spin model 
are assumed to take those expected for the 3D classical 
Heisenberg model. The phase-transition point, Ac, is 
then determined by the formula XAF ~ (Ac — A)"''' and 
A ~ (Ac — A)'' with the critical exponents 7 = 1.4 and 
u = 0.71.17) 

§4. Quantum Phase Transitions 

We now discuss quantum phase transitions when the 
competing exchange couplings are varied. There are 
three parameters in the perturbation term. A, a and k. 
We first study the model given by Fig. 1(a), which con- 
tinuously connects the 2D Shastry-Sutherland model and 
the ID orthogonal-dimer chain, and then move to the 
model (b), which includes the Heisenberg model on the 
1/5 depleted square lattice. 
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4-1 Adiabatic evolution from the orthogonal- dimer 

chain 

Let us start with the system shown in Fig. 1(a). Note 
that this spin system is reduced to the orthogonal-dimer 
chain for /c = 0, which was aheady studied in detaih 
when the ratio a of the competing exchange couphngs is 
varied, the first-order quantum phase transition between 
the spin-gap phases with the dimer and the plaqucttc 
structures occurs at ad = 0.81900.^^'^*^ The merit to 
study this system is that we can check how the known 
spin-gap states in the chain system adiabatically evolve 
when the system approaches the 2D Shastry-Sutherland 
model. 

We first discuss the effect of interchain coupling k 
on the first-order phase transition between two spin-gap 
phases by computing the ground state energy Eg. Note 
that the energy for the dimer phase is unchanged even on 
the introduction of k due to the orthogonal-dimer struc- 
ture. For the plaqucttc phase, by choosing the isolated 
plaquettes shown by shaded squares in Fig. 1(a) as an 
initial configuration, we compute the ground state en- 
ergy Eg up to the sixth order in A with a and k being 
fixed. By applying the first-order inhomogeneous differ- 
ential method to the obtained series, we then estimate 
the ground state energy for the plaqucttc phase by set- 
ting A = 1. The energy estimated for k = 0.1, 0.5 and 
0.9 is shown in Fig. 2. It is seen that the introduction 
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Fig. 2. Ground state energy per site as a function of a = J(/J 
(for A = 1). The flat line {Eg/JN = -3/8) is the energy for the 
exact dimer state. The dotted, dashed and solid lines denote the 
energy at A = 1 obtained by the plaquette expansion for k = 0.1, 
0.5 and 0.9, respectively. 

of k prefers the plaquette phase, shifting the first-order 
transition point ad to the left in Fig. 2. From these 
analyses, we determine the phase boundary between the 
dimer phase and the plaquette phase in the 2D param- 
eter space, which is shown in Fig. 3. We have checked 
that the present estimate of the critical point for fc = 
reproduces the above-mentioned value of ad fairly well. 
It is seen that the phase boundary is not so sensitive to 
k for small k, because the energy for the plaquette phase 
decreases qtiadratically in k. 

It is to be noted that these two competing spin- 
gap phases are continuously connected to those in the 
Shastry-Sutherland model (fc = 1), with the nature of the 
first-order transition being unchanged. To check the va- 
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Fig. 3. Phase diagram for the 2D frustrated spin system shown in 
Fig. 1(a). The phases I, II and III represent the dimer, plaquette 
and magnetically ordered phases, respectively. The dotted line 
corresponds to the 2D Shastry-Sutherland model ajid the dot- 
dashed lines to the asymptotic phase boundary for large a. The 
cross indicates the parameters deduced for SrCu2 (603)2.^®' 



lidity of the obtained results, we have also performed the 
plaquette expansion starting from another initial config- 
uration, where an isolated plaquette is formed by four 
spins shown as 1, 2, 3 and 4 in Fig. 1(a). From the 
ground state energy computed up to the seventh order, 
we confirm that the phase boundaries obtained by two 
distinct ways arc indeed consistent within our accuracy. 

To complete our discussions on this model, it is nec- 
essary to discuss whether the system is driven to the 
magnetically ordered phase with the increase of the in- 
terchain coupling. To this end, we calculate the stag- 
gered susceptibility xaf and the spin gap A by means 
of the plaqucttc expansion up to the fourth and the fifth 
order in A, respectively. We show the results obtained 
for k = 0.9 and 0.8 in Fig. 4. It is seen that two lines de- 
termined from different quantities provide consistent val- 
ues even around A = 1, implying that the phase bound- 
ary may be obtained rather accurately for the original 
model (A = 1) although our calculation is restricted to 
the lower-order series expansion. From the above analy- 
ses, we end up with the phase boundary shown in Fig. 3, 
which separates the plaqucttc phase II and the magneti- 
cally ordered phase III. Note that for large a, the bound- 
ary between II and III approaches the correct asymptotic 
line, which was estimated as fc,, = 0.56.^"' This also sup- 
ports that the present calculation gives the reliable re- 
sults for the phase diagram. 

In the obtained phase diagram shown in Fig. 3, we 
remark again that two competing spin-gap states in 
the orthogonal dimer chain^^'^*' are continuously con- 
nected to those in Shastry-Sutherland model. In partic- 
ular, a RVB-like intermediate spin-gap phase^^ proposed 
for the Shastry-Sutherland model naturally emerges via 
the adiabatic continuation of the plaquette phase in the 
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Fig. 4. Second-order phase transition between the plaquette 
phase and the antiferromagnetically ordered phase for the 2D 
spin system given by Fig. 1(a) (for k = 0.9 and 0.8). The solid 
(dashed) line indicates the phase boundary obtained by biased 
[3/2] ([2/2]) Fade approximants for the spin gap (the staggered 
susceptibility) . 



orthogonal-dimer chain. 

4-2 Adiabatic evolution from the 1/5-depleted square 
lattice 

We next deal with the spin system with a shghtly dif- 
ferent structure shown in Fig. 1 (b). Recall that when 
k = 1 [k = 0), this model with A = 1 is reduced to 
the Shastry-Sutherland model (1/5-depleted square lat- 
tice model). We here study how the spin-gap phase 
in the orthogonal-dimer model for SrCu2 (603)2 is re- 
lated to that in the 1 /5-depleted square-lattice model for 
CaV409, when the frustrating interaction k is continu- 
ously varied. We carry out the dimer (plaquctte) expan- 
sion for the ground-state energy, the staggered suscepti- 
bility and the spin gap. By performing similar asymp- 
totic analyses mentioned above, the phase diagram is 
determined, which is shown in Fig. 5. The detail of 
calculation will be described below for each case. 

Let us start our discussions with the case of fc = 0, i.e. 
the 1/5-depleted square lattice model, which was already 
studied in connection with CaV409.^^-' In this 
the ratio of the exchange interactions a is changed, the 
dimer phase I first undergoes the second-order phase 
transition to the antiferromagnetically ordered phase IV 
at the critical point a^2 = 0.590(5), and this ordered 
phase further shows the second-order phase transition 
to the plaquette phase II at = 1-09(1). The criti- 
cal points obtained by the present method agree fairly 
well with those of quantum Monte Carlo simulations. 
Note that the antiferromagnetically ordered phase on 
1 / 5-depleted square lattice labeled by IV is different from 
that on the square lattice indicated by III. 

As the coupling k is increased, the spin-gap phases 
become more stable while the magnetic correlation for 
the phase IV is suppresscid, whic;h is mainly due to the 
frustrating nature of the interaction k. It is seen that 
the competition among three phases persists up to the 
multicritical point (a, ka) ~ (0.72, 0.32) beyond which 
the magnetically ordered phase IV disappears. Thus for 
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Fig. 5. Phase diagram for the 2D frustrated spin system shown 
in Fig. 1(b). The phase I and II represent the dimer and the pla- 
quette phase. The phase III (IV) represents the antiferromagnet- 
ically ordered phase on the square lattice (1/5-depleted square 
lattice). The dotted line indicates the 2D Shastry-Sutherland 
model and the dot-dashed lines the asymptotic phase bound- 
ary for large a. The crosses + and x indicate the location of 
SrCu2(B03)2i^) and CaV409,9) respectively. 



larger k the dimer state directly undergoes the first-order 
phase transition to the plaquette spin-gap phase. Note 
that the dimer phase here does not have the exact eigen- 
state except for the orthogonal-dimer case fc = 1, which 
is contrasted to the model (a). We thus have estimated 
the ground state energy both of the dimer and the pla- 
quette phases in this case;. For reference, we show in Fig. 
6 the ground-state energy calculated for k = 0.9, where 
the dimer (plaquette) expansion has been first performed 
up to the eighth (seventh) order in A, and then the first- 
order inhomogeneous differential method has been used 
to obtain the values at A = 1. To confirm the valid- 
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Fig. 6. Ground state energy per site as a function of o = J( / J 
(for A = 1 and k = 0.9) for the model shown in Fig. 1(b). The 
solid (dashed) line indicates the energy for the dimer (plaquette) 
phase. We also show the results obtained by the exact diagonal- 
ization for AT = 24 sites with a periodic boundary condition as 
open circles. 



ity of our series expansion, we have also performed the 
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exact diagonalization for the system with N = 24 sites 

with a periodic boundary condition, which is found to be 
consistent with the series-expansion results. The phase 
boundary thus obtained between two spin-gap phases is 
drawn in Fig. 5. We note that our series expansion turns 
out to be not so accurate around the tricritical point, so 
that the phase boundary around there is shown as the 
dashed Une. 

When the coupUng k is further increased, the plaquctte 
phase shows the instabihty to an antiferromagnctically 
ordered phase III on the square lattice. To study this 
instability from II to III, we have performed the plaquc- 
tte expansion for the staggered susceptibility (spin gap) 
up to the fourth (fifth) order, and dctcrmincxi the phase 
boundary using the biased Fade approximants. which is 
shown in Fig. 5. 

Consequently, we arrive at the phase diagram shown 
in Fig. 5, in which the parameters deduced for the spin- 
gap compounds CaV409 and SrCu2 (603)2 arc also in- 
dicated. It is clearly seen how the 2D Heisenberg model 
on the 1/5-depleted square lattice is connected to the 
Shastry-Sutherland model. It is particularly interesting 
to notice that a RVB-type spin gap state proposed for 
the Shastry-Sutherland model is adiabatically connected 
to a plaquette singlet state for CaV409. 
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Fig. 7. Spin-triplet excitation energy A/ J for (a) Shastry- 
Sutherland model and (b) orthogonal-dinicr chain as a function 
of tt(= J[/J). In (a), a triplet excitation spectrum^'*^ over the 
dimer state is plotted for a < ad, while two lowest triplet ex- 
citations are shown for the RVB-type phase ad < a < a.ci- In 
(b), similar plot is done for the dimer phase (a. < ac) and also 
for the plaquette phase (ac < a).^^) The meaning of the solid, 
dashed and dotted lines is mentioned in the text. 



4-3 Excitation spectrum 

We have seen that the spin-gap phases in the Shastry- 
Sutherland model are continuously connected to the 
spin-gap phases known for two different models. Ac- 
cording to this relationship, it may be expected that the 
Shastry-Sutherland model shares some of its character- 
istic features with the latter models. As an example, we 
here show that the spin-triplet excitation spectrum for 
the RVB-type phase in the Shastry-Sutherland model in- 
deed exhibits a similar behavior expected for the plaque- 
tte phase in the orthogonal-dimer chain. In Fig. 7(a), 
we show the spin-triplet excitation spectrum calculated 
for the 2D Shastry-Sutherland model. In contrast to the 
dimer phase, where the lowest triplet excitation is an 
almost-localized triplet over the dimer ground state,'') 
we find that the lowest-triplet excitation in the RVB-type 
phase shows a more complicated feature; two triplet exci- 
tations show the level-crossing around value of a' ~ 0.76. 
It is not easy to clearly see the origin of this level-crossing 
because our calculation for the 2D model is restricted to 
lower-order perturbation. Fortunately, we can exploit 
the idea of the adiabatic continuation by continuously 
changing the system to the orthogonal-dimer chain. 

In the orthogonal-dimer chain, a similar level-crossing 
was recently observed,^^) as shown in Fig. 7(b). In this 
case, the nature of two different excitations in the pla- 
quette phase is well understood. Namely, in the lowest 
excitation denoted by the dotted line in the region of 
a > a'(~ 0.87) is a triplet excitation that simply breaks 
a plaquette singlet. On the other hand in the region 
ttci < a < a', the lowest-energy excited state shown 
by the dashed line is described by a four-fold degenerate 
level without dispersion, in which two diagonal spins 
on a square plaquette becomes completely free, giving 
rise to the four-fold degeneracy. The emergence of this 



peculiar excitation reflects the fact that the system is lo- 
cated closely to the first-order phase transition point. 
In other words, this level crossing is due to the strong 
frustration effect around the transition point. When the 
interchain coupling k is introduced, the above two-types 
of excitations are continuously changed to those in the 
Shastry-Sutherland model, as shown in Fig. 7(a). In 
particular, four-fold degenerate states in the chain sys- 
tem split into singlet and triplet states, and the latter 
gives the lowest-excited state shown by the dashed line 
in the Shastry-Sutherland model in Fig. 7(a). In this 
way, the characteristic level-crossing found in the ID sys- 
tem persists even for the 2D system, from which we can 
say that the unusual level-crossing found in the Shastry- 
Sutherland model reflects the strong frustration around 
the first-order phase transition point. 

§5. Summary 

We have studied the ground-state phase diagram for 
a frustrated quantum spin model on the square lattice, 
which systematically describes the systems for the spin- 
gap compounds such as SrCu2(B03)2 and CaV409. By 
calculating the ground state energy, the staggered sus- 
ceptibility and the spin gap by means of the series ex- 
pansion method, we have determined the phase diagram, 
and clarified the nature of the associated quantum phase 
transitions. One of the main purpose of the present study 
is to uncover the origin of the intermediate spin-gap state 
proposed for the Shastry-Sutherland model, which has a 
spatially extended nature. Starting from two types of 
the well-studied models, we have observed how this in- 
termediate phase emerges in the phase diagram. Namely, 
with the use of the simple orthogonal-dimer chain model, 
it has been checked that the first-order phase transition 
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between the dimer and RVB-type phases in the Shastry- 

Sutherland model has the same origin as in the transi- 
tion between the dimer and plaqucttc phases known for 
the chain system. Further analysis based on the second 
model has allowed us to show that the RVB-type inter- 
mediate state is adiabatically connected to the plaquctte 
state known for the 1/5-depletcd square lattice model, 
which describes the essential properties of CaV409. We 
have also shown that the level-crossing of the triplet ex- 
citations in the Shastry-Sutherland model is understood 
well in terms of their counterparts in the ID orthogonal- 
dimer chain model. 

These systematic investigations based on the adia- 
batic continuation support that the intermediate spin- 
gap phase for the Shastry-Sutherland model may be sta- 
bilized by geometrical frustration induced by the compet- 
ing interactions. However, we wish to mention that there 
still remain unresolved problems for the phase diagram 
of the; Shastry-Sutherland model. For example, Weihong 
et al. have shown by means of the Ising expansion^**) 
that the system has the finite staggered magnetization 
down to {J' /J) — 0.71, which may contradict the fact 
that the spin gap phase is extended up to (J'/ J) ~ 0.86. 
Although we have checked that the ground state energy 
of the present spin-gap phase for [J' / J) < 0.86 is in- 
deed lower than that of Weihong et al., this apparent 
inconsistency for the magnetization should be resolved 
in the future study. Also, Knetter et al. have recently 
claimed that the instability of the two-magnon excitation 
in the dimer phase may occur at {J'/ J) = 0.630(5).^^ Al- 
though it may not be clear whether this indeed triggers 
a first-order phase transition to some new phase, it is 
an important open problem to answer the above points 
for the phase diagram of the Shastry-Sutherland model. 
This is now under consideration. 
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